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The application of Wilson's Numerical Renormalization Group (NRG) method to dissipative 
quantum impurity models, in particular the sub-ohmic spin-boson model, has led to conclusions 
regarding the quantum critical behavior which are in disagreement with those from other methods 
and which are by now recognized as erroneous. The errors of NRG remained initially undetected 
because NRG delivered an internally consistent set of critical exponents satisfying hyperscaling. 
Here we discuss how the conspiracy of two errors - the Hilbert-space truncation error and the 
mass-flow error - could lead to this consistent set of exponents. Remarkably, both errors, albeit of 
different origin, force the system to obey naive scaling laws even when the physical model violates 
naive scaling. In particular, we show that a combination of the Hilbert-space truncation and mass- 
flow errors induce an artificial non-analytic term in the Landau expansion of the free energy which 
dominates the critical behavior for bath exponents s < 1/2. 
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I. INTRODUCTION 

The Numerical Renormalization Group method^' 2 - 
originally developed^ by Wilson for the Kondo model and 
subsequently applied to a variety of impurity models with 
fermionic baths, was generalized in Refs. |3||j to the case 
of bosonic baths with an eye towards dissipative impu- 
rity models. In particular, it has been applied to the 
spin-boson model, 5 



where a z = ±1 are the local impurity states with bias e, 
f2 is the tunneling rate, and the LOi > and Ai are the 
frequencies and coupling constants of the bath oscillators. 
The bath is completely specified by its propagator at the 
"impurity" location 



r( w ) 



E 



A 



+ i0+ 



(2) 



whose spectral density J(ui) = 
parameterized as 



Tmr(w) is commonly 



J{lo) — 2tt a <jj\ s lu s , < u) < uj c 



s > -1 



(3) 



where the dimensionless parameter a characterizes the 
dissipation strength, and uj c is a cutoff energy. The 
value s = 1 represents the case of ohmic dissipation. 
For < s < 1, the spin-boson model ([1]) displays a 
quantum phase transition (QPT) between a delocalized 
and a localized phase, reached for small and large a, 
respectively i^r— 

NRG has enabled detailed investigations of this QPT. 
While statistical-mechanics argument a 8 i 10 ' 16 suggest that 
this transition is in the same universality class as the 
thermal phase transition of the one-dimensional (Id) 
Ising model with l/r 1+s long-range interactions, initial 



NRG results 7 were in disagreement with this quantum- 
to-classical correspondence (QCC). In particular, NRG 
delivered s-dependent critical exponents which obeyed 
hyperscaling for all < s < 1. Those exponents agreed 
with the ones of the Ising model for s > 1/2, but dis- 
agreed for s < 1/2 where the Ising model is above its 
upper-critical dimension and displays mean-field behav- 
ior without hyperscaling . 17 ' 18 This apparent violation of 
QCC prompted further numerical investigations. Quan- 
tum Monte Carlo (QMC)^ and exact-diagonalizationii 
studies led to the opposite conclusion, namely that the 
critical behavior of the spin-boson model for s < 1/2 is 
classical and of mean-field type. 

Meanwhile, two different errors of the NRG method 
have been identified, 19 which can be held responsible for 
incorrect results regarding critical properties. These two 
errors are the Hilbert-space truncation error, arising from 
the fact that the infinite Hilbert space of each bath har- 
monic oscillator needs to be truncated to Nb local states, 
and the mass-flow error, which is rooted in the itera- 
tive diagonalization procedure of NRG and leads to an 
effective shift of model parameters along the NRG run. 
The mass-flow error has been discussed in some detail in 
Ref. [13 where it has been shown to give rise to an incor- 
rect exponent x, describing the thermal divergence of the 
susceptibility at criticality. Ref. [2(] also implemented a 
recipe to correct the mass-flow error which then allowed 
to recover the mean-field value x = 1/2 for s < 1/2. The 
problem of Hilbert-space truncation has been solved in 
Ref. [l3|, by employing an optimized boson basis within 
a variational matrix-product-state approach on the Wil- 
son chain. Applied to the spin-boson model with s < 1/2, 
this method not only delivered the expected mean-field 
values for the order-parameter exponents f3 and S if the 
variational bosonic Hilbert space is taken to be (effec- 
tively) infinite, but also showed that artificially restrict- 
ing the Hilbert space causes a crossover at small energies 
to the erroneous non-mean-field behavior of NRG. 
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The shortcomings of NRG have independently become 
clear in studies of dissipative anharmonic-oscillator mod- 
els: These models are in one-to-one correspondence with 
a Id Ising </> 4 theory with l/r 1+s interactions and must 
therefore display a mean-field transition for s < 1/2. 
Nevertheless, standard NRG delivered^ the same non- 
mean-field exponents as for the spin-boson model. 

Together, these developments show that the initial 
NRG results^ were incorrect and that the QCC is fulfilled 
for the sub-ohmic spin-boson model. However, it has re- 
mained a puzzle why and how the NRG could produce 
an internally consistent set of exponents which moreover 
obeyed hyperscaling. This question is of obvious rele- 
vance for further applications of NRG. It is the purpose 
of the present paper to close this gap. We discuss Hilbert- 
space truncation in some detail, showing that the com- 
bination of Hilbert-space truncation and mass flows acts 
like a non-analytic term in the Landau expansion of the 
free energy. At criticality, this in turn forces the system 
to effectively remain locked at the upper-critical dimen- 
sion for all s < 1/2. As a result, a set of "trivial" ex- 
ponents with hyperscaling properties is obtained. While 
Hilbert-space truncation and mass flow may be thought 
of as independent sources of error, their effects are in- 
tertwined: As we argue below, the non-analytic term in 
the Landau expansion can be understood in terms of a 
truncated mass flow. 

We note that recent papers^— insist on the break- 
down of QCC for the spin-boson model. However, the 
concrete numerical results obtained in Refs. [2lll23l are 
fully consistent with the scenario proposed here, as will 
be discussed below. Moreover, the objections raised 
against QCC in Ref. do not apply. Hence, we believe 
the interpretation given in Refs. I2ll423l is incorrect. 



A. Outline 

The remainder of the paper is organized as follows: 
Sec. Illl summarizes the spin-boson model, together with 
the standard derivation of its QCC via a Feynman path 
integral. It then reviews the properties of the resulting 
4> 4 theory - whose language will be heavily used in the 
paper - in particular the mean-field regime of s < 1/2 
where Landau theory applies. Last but not least, Sec. |TT] 
proposes that supplementing the field theory by a cer- 
tain non-analytic term yields exactly the non-mean-ficld 
exponents seen in NRG - this non-analytic term will 
derived and discussed in the remainder of the paper. 
Sec. Mil dives into technical details of NRG, by repeat- 
ing and extending the discussion of the mass-flow error 
as given in Ref. [20 - this includes sub-leading correc- 
tions as addressed in Ref. [H. Finally, Sec. [TV] is devoted 
to the Hilbert-space truncation error and the truncated 
mass flow. It will highlight the crossover scales induced 
by the truncation and their influence on physical observ- 
ables. A general discussion of NRG and its applicability 
to mean-field critical phenomena closes the paper. 



II. SPIN-BOSON MODEL AND FIELD THEORY 

On general grounds, one expects that the quantum 
phase transition in the spin-boson model can be described 
by a field theory for an Ising (i.e. real scalar) order pa- 
rameter 4>{t), where (0) corresponds to the local magne- 
tization (<7 Z ). 



A. Derivation 

Such a field theory can be derived from the Feynman 
path-integral representation^ of the partition function 
of %sb- One starts with a Trotter decomposition of the 
imaginary time axis into N intervals of length At = fi/N. 
The operators exp(— 'HsbAt) are then evaluated by in- 
serting identity operators in terms of the eigenstates of 
o z , with eigenvalues S = ±1, and the bath-oscillator po- 
sitions, qi, for each time slice k. (Utilizing coherent states 
is not necessary, as Hsb contains distinguishable objects 
onlyj2£) One arrives at: 
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where m, are the masses of the bath oscillators, and the 
limit At — > has to be taken at the end as usuali^i 
Note that the action of a x in 'Hsb has been evaluated 
directly. The oscillator bath is best treated in Fourier 
(i.e. frequency) space. Defining 
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for integer wavenumbers n = 0, . . . , N — 1 and using the 
explicit form of the spectral density ([3]) , one finds for the 
bath term the exact result: 



(^Ar) 1 - 8 
sin(7rs/2) 



2 sin — 

N 



(6) 



The integrals in partition function involving the bath- 
oscillator coordinates are Gaussian, such that the bath 
oscillators can be integrated out exactly^ After trans- 
forming back to real space (i.e. time) the full action takes 
the form 

N 1 

S = -^[-lncoth(ftA-r)]S fc Sfc +1 (7) 

k=l 

+ l^2K(k-k')(S k -S k ,) 2 . 
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The bosonic bath thus generates a long-ranged self- 
interaction in imaginary time, with the long-distance be- 
havior 

K(k-k') = oT(l + s)( Uc ^-' {k _l )1+a , (8) 

valid for 1 <C \k— k'\ <C N. The representation is real 
and free of a Berry-phase term^ 

At this stage, one may decide to interpret the Sk (for 
fixed finite At) as discrete spins, which form an Ising 
chain with long-range interactions. Its thermodynamic 
limit is obtained as /3 — ¥ oo, and one can show that the 
long-distance properties are independent of At. (This is 
somewhat non-trivial, as the nearest-neighbor coupling 
in ||7j! diverges as At — > 0, while the long-range coupling 
vanishes in the same limit.) Speculations that the limits 
At — > and (3 — > oo might not commute^ have turned 
out to be unwarranted, as has been shown by explicit 
numerical computations^ The phase transition in this 
Ising model itself is described by a continuum 4> 4 theory 
with long-range interactions^ 

Alternatively, one may depart from Eq. (J3) and relax 
the hard-length constraint on the Ising variables. Assum- 
ing a slowly varying Ising order-parameter field <f> then al- 
lows to obtain a continuum theory with finite coefficients, 
which is - of course - the same <p A theory as obtained for 
the discrete Ising model. 

This 4 theory reads 

S = J ^(mo + M*)W<w)| 2 + J dr K(t)+#M] 

(9) 

the continuous frequency integral has to be replaced by a 
Matsubara sum for T > 0. The field <j> has been re-scaled 
such that the prefactor of the |w| s term equals unity, and 
e oc e represents the field conjugate to the order parame- 
ter. In Eq. © an increase of the (bare) mass parameter 
mo corresponds to a decrease in the dissipation strength 
a, with the critical point located at mo — m c or a = a c . 

By universality arguments, the same field theory de- 
scribes a dissipative anharmonic oscillator 

"Hdho = Cla i a+^(a + at) + Un a (n a -1) (10) 
+ 2S A *(a + a t )(&i + ^ t )+S w * 6 * 6 *' 

i i 

where > is the bare "impurity" oscillator frequency, 
n a = a^a, and U parameterizes the anharmonicity. For 
U = oo, 'Hdho is equivalent to Hsb, because all impurity 
states except n a =0,1 are projected out. In the opposite 
limit of small U the interaction term of 'Hdho (HHJ) is 
expected to follow the same RG flow as does u in S ©• 



B. Critical exponents 

The field theory ([9]) can be utilized to calculate critical 
exponents, with the following definitions: 

<f)(a > a c ,T = 0,e = 0) oc {a-a c f, 

x(a<a c ,T = 0) oc (a c -a)" 7 , 

4>(a = a c ,T = 0) oc lei 1 / 5 , 

X (a = a c ,T) oc T~ x , 

x"(a = a c ,T = 0,w) oc M^sgnH, 

where <f> = ((f>) and x — dcj>/de. The last equation de- 
scribes the dynamical scaling of x- 

The 4> 4 theory ([§]) has been analyzed in detail in 
Ref. Il7l . Power counting yields the scaling dimensions 
at criticality: 

din#(T)] = (l-s)/2, (12) 
dim[u] = 1 — 4dim[0(T)] = 2s — 1 , 

i.e., the system is above (below) its upper-critical dimen- 
sion for s < 1/2 (s > 1/2). Consequently, the transition 
is controlled by a Gaussian fixed point for s < 1/2, with 
critical exponents 

Pmf = 1/2 , 7MF = 1 , *MF = 3 , V M F = 1/s , (13) 

where the s dependence of the correlation-length expo- 
nent v reflects the influence of the long-range interac- 
tions. Moreover, the exponents x and y, related to the 
finite-size-scaling and anomalous decay exponents of the 
classical Ising model, are different: 

IMF = 1/2 , 2/mf = s . (14) 

In contrast, for s > 1/2 the exponents take non-trivial 
values. Hyperscaling allows to deduce a few exact results, 
given that y = s is known to be exact i 18 i 27 

1 + s 

•^hyp — yhyp = s j ^hyp = Z ■ (1^) 



C. Temperature-dependent mass and susceptibility 

As a prerequisite for the following considerations, we 
now explicitly derive the temperature dependence of the 
(renormalized) order-parameter mass at criticality, which 
is related to the susceptibility by 1/%(T) = m(T). 

Below the upper-critical dimension, i.e., for s > 1/2, 
naive scaling applies. Given the zero-temperature form 
of the critical propagator l/G^ oc uj s , the interaction- 
induced mass at criticality has to scale as m(T) oc T s , 
which yields x = s as stated above. 

In contrast, in the mean-field regime of s < 1/2, m(T) 
can be obtained in lowest-order perturbation theory in 
u. The tadpole diagram at criticality evaluates to 

= +6 j (16) 

7 <* 
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where ujj — 2-kjT are Matsubara frequencies and bojj, 
representing the short-range piece of the interaction, 
has been included for convergence. The temperature- 
dependent piece, £i(T) — Ei(0) = m(T), can be obtained 
by simple scaling analysis. Assuming m(T) = aT x we 
find 



aT x = uT 1 ' 
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(17) 

For x > s and small T the right-hand side is dominated 
by the j — piece of the sum and hence approaches a 
constant as T —> 0. Demanding the temperature power 
laws on both sides of Eq. ([T7| to match, we obtain x — 
1/2 as expected^! (which indeed fulfills x > s). 

The calculation shows that it is mandatory to include 
the mass self-consistently into the internal propagator of 
the tadpole diagram; this is different from the situation 
in c/) 4 theories with short-ranged interactions where the 
tadpole diagram with bare propagator would be sufficient 
to determine the shift exponent. 



D. Field theory with NRG errors 

Our central result is that the errors of NRG induce ex- 
tra non-analytic mass terms in the order-parameter field 
theory: 



5 orr = / dr0 2 (r) IvTo + wM] 2 ^ 1 -^ 
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with v and w depending on parameters of the NRG algo- 
rithm. The first term is a consequence of the mass flow, 
as discussed in Ref. [20) and summarized in Sec. lIIII below. 
It is responsible for the incorrect result x = s found us- 
ing NRG for s < 1/2. The second term originates from 
Hilbert-space truncation, with the prefactor w scaling 



with the number Nh of bosons as 



n: 



s/(l- 



It will 



be derived in detail in Sec. IIVI 

Let us discuss the effect of the second term on critical 
exponents. Approximating 4>(t) as static in iS + iS err and 
setting T = leads to a Landau free energy 



F(<j)) = e(f> + m4> 2 + w\4>\ 2/(1 - s) + ucj) 4 



(19) 



It is important to note that fluctuation corrections to 
this mean-field picture arise from the quartic u term only 
(as w is not an interaction), and hence will influence the 
critical behavior only for s > 1/2, where moreover 2/(1 — 
s) > 4 and hence the w term is subleading. In contrast, 
for s < 1/2 simply minimizing F(<j>) (|19[) is sufficient to 
calculate the zero-temperature exponents. For w = this 
of course results in /3mf = 1/2 and <5mf = 3. However, a 
non-zero w term dominates over the quartic u term for 
s < 1/2, and minimizing F(<p) then gives 
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FIG. 1: (Color online) Structure of the NRG Hamiltonian, 
with the bath represented by a semi-infinite Wilson chain, 
with bath operators b' n . The boxes indicate the iterative di- 
agonalization scheme. 



which are exactly the exponents observed^ in NRG! The 
equation of state derived from (|19p is also consistent with 
the comprehensive scaling analysis of NRG results at dif- 
ferent Nb performed in Ref. u&i 

In the remainder of the paper, we discuss the errors 
of the bosonic NRG and the consequences of Eq. (|T8| in 
detail. 



III. NRG AND MASS FLOW 

The mass-flow problem of the bosonic NRG has been 
investigated in Ref. ;2C| Here we shall summarize the dis- 
cussion and also provide an analysis of certain subleading 
corrections. 

Within the NRG algorithm, the bath is represented by 
a semi-infinite ("Wilson") chain, Fig. [TJ such that the 
impurity is coupled to the first site of the chain only, and 
the local density of states at this first site is a discrete ap- 
proximation to the bath density of states*^ Due to the 
logarithmic discretization, the site energies e n and hop- 
ping matrix elements t n decay exponentially along the 
chain according to w c A~ n+1 , where A is the discretiza- 
tion parameter. We define a Hamiltonian % n for the 
impurity plus the first n sites of the Wilson chain, and 
we denote by (w) the corresponding propagator at the 
impurity site of this n-site bath. Then, Hoc with bath 
roo(w) is the discretized version of the original problem. 

NRG proceeds by iteratively diagonalizing the Hamil- 
tonian: First, Hi is diagonalized and the lowest N s eigen- 
states are kept. Then, the next bath site is added to form 
%2j the new system is diagonalized, and again the lowest 
N s eigenstates are kept (which are approximations to the 
lowest states of T^)- As the characteristic energy scale of 
the low-lying part of the eigenvalue spectrum decreases 
by a factor of A in each step, this process is repeated un- 
til the desired lowest energy is reached. Thermodynamic 
observables at a temperature T n = w c A~' i+1 //3 are typ- 
ically calculated via a thermal average taken from the 
eigenstates at NRG step n. Here, j3 is a parameter of 
order unity which is often chosen as j3 = 1. 



A. Mass-flow error 

The iterative diagonalization procedure implies that, 
at NRG step n, the chain sites n + 1, n + 2, ... have 
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not yet been taken into account, i.e., the effect of those 
sites does not enter thermodynamic observables at tem- 
perature T„. The "missing" sites imply a missing con- 
tribution to the real part of the bath propagator. As 
detailed in Ref. |2Cj, the zero- frequency value of this miss- 
ing real part, Re(r oo — T n )(uj — 0), scales as (f3T n ) s for 
a power-law bath spectrum, Eq. ([3]). In the field-theory 
language, the missing real part implies a temperature- 
dependent variation of the order-parameter mass. This 
is most clearly seen for the dissipative oscillator model 
in Eq. (|10|) where Rer(u; = 0) directly renormalizes the 
oscillator frequency, fl — > + Rer(u; = 0) at U = (note 
that ReT(w = 0) < 0). Similarly, for the spin-boson 
model it is clear that the missing real part effectively 
leads to a temperature-dependent variation of model pa- 
rameters. Provided that the spin-boson model renormal- 
izes to the same (j> 4 field theory as the dissipative oscil- 
lator, as shown in Sec. |4j it is clear that this variation is 
in the (j> mass. 

Thus, the mass- flow error generates an artificial mass 
contribution m a = vT s - this is the first term in <S crr (|18|) . 
The prefactor v scales as /3 s , but also depends on the 
model and parameter regime under consideration (note 
that the mapping of the microscopic H to the field theory 
involves coarse graining). Moreover, interaction effects 
will generate subleading contributions to m a , see below. 

B. Temperature dependence of \ m the presence of 
mass flow 

As derived in Sec. Ill CI the irrelevant interaction u 
leads to an order-parameter mass oc T 1 / 2 along the flow 
towards the Gaussian fixed point for s < 1/2. Conse- 
quently, the physical susceptibility follows x T~ x / 2 
(Ref. \17\ ). However, the artificial mass m a oc T s domi- 
nates the physical mass at low T, leading to the unphys- 
ical result l/\ oc T s - which happens to coincide with 
the physical result for an interacting fixed point with hy- 
perscaling. 

Considering the results in Ref. [23|, it is worth eval- 
uating the subleading correction to the mass-flow- 
determined m a = vT s . This piece, rri2(T), is now aris- 
ing from the quartic interaction and can be evaluated in 
lowest-order perturbation theory. In analogy to Eq. (TT5)) 
we need to evaluate Si(T) — Si(0) = m2(T), with 

Ei = uT V : — (21) 

^ \u n \ s + bcol + vT s + m 2 {T) y ' 

where the propagator now contains both the artificial and 
the interaction-generated mass contribution. Assuming 
m^T) = V2T X2 with X2 > s, the term in the inter- 
nal propagator is subleading. Consequently, the equation 
determining m,2 reads 

V2T 2 = UT ^ (| 2^.7 I s + v)T s ~ J \2ttj\ s T^ (22) 



and one immediately arrives at x 2 — 1 — s. Hence, the 
presence of the artifical mass modifies the interaction- 
generated mass.— Indeed, the NRG results for x(T) for 
both the spin-boson model 23 and the dissipative anhar- 
monic oscillator™ are consistent with X2 = 1 — s and 
thus in line with the above calculation. The calculation 
falsifies the claim of Ref. [23| that X2 > 1/2 would be 
incompatible with the mass-flow scenario. 



C. Mass-flow correction 

Ref. HO proposed an empirical algorithm to cure the 
mass-flow problem. This is based on compensating 
the NRG-induced mass flow by including an explicit 
temperature-dependent term, with coefficient kT s , in the 
Hamiltonian which shifts the critical point similar to the 
effect of ReT. For the spin-boson model, the natural 
choice is to add a temperature-dependent piece to the 
tunneling matrix element fi. In this procedure, K re- 
mains a free parameter, and Ref. [2^ proposed to adjust 
its value at criticality such that the unphysical T s piece 
in 1/x is removed. This algorithm, dubbed NRG*, was 
shown to work well and to yield the expected mean-field 
behavior 1/x oc T 1 / 2 at low temperature for s < 1/2 
whereas the T s behavior was stable for s > 1/2. There- 
fore, the QCC for the exponent x has been confirmed by 
NRG*. 

Ref. [23] recently questioned this result, based on the ob- 
servation that l/x(T) from NRG*, which follows T 1 / 2 at 
low T, displays an unusual downward-deviation from this 
law at higher T (see Figs. 7b and 8b of Ref. HO). Ref. [H 
claimed that this would reflect an underlying temper- 
ature power law with exponent smaller than 1/2, then 
supposedly inconsistent with mean-field behavior. Here 
we argue instead that the downward-deviation can be 
straightforwardly understood in terms of the mass flow 
and its compensation. As explained in Ref. the mass- 
flow correction term with fixed k only compensates for 
the mass flow near a specific fixed point, and different 
fixed points would require different k. In other words, 
if the proper kq for the low-temperature fixed point has 
been chosen, then this will compensate the artificial mass 
only to leading order, but a subleading term will be left 
over. For the spin-boson model, one can easily under- 
stand the sign of the subleading term by noting that 
NRG* produces an overcompensation of the mass flow 
at elevated temperatures or energies: At high energies, 
ReT influences the mass very little (note that it cou- 
ples to cr 2 in Eq. (|TJ) in the short-time limit). Then, 
at high energies, only the mass-flow correction is active, 
which leads to a downward correction of the mass. It 
is this downward correction that causes the downturn 
in l/x(T). (We also note that the downturn does not 
correspond to a new power-law regime, contrary to the 
assertion of Ref. |23j.) We conclude that the downturn in 
l/x(T) is not physical, but instead related to the short- 
comings of the simple mass-flow correction in NRG* . The 
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same downturn is observed when applying NRG* to the 
dissipative anharmonic oscillator with large U . 



IV. NRG AND HILBERT-SPACE TRUNCATION 

In this section, we discuss how the truncation of the 
bosonic Hilbert space at each site of the Wilson chain, 
i.e., using only the lowest Nb occupation- number eigen- 
states of the oscillator as basis, influences the evaluation 
of observables in the spin-boson model. In particular, we 
will explain how the non-analytic w term in the effective 
field theory (fTHJ) emerges. 

Hilbert-space truncation is bound to become impor- 
tant if <r z develops a finite expectation value, either inside 
the ordered phase or in an applied field e. For a finite 
(<r z ), the bath oscillators will have a finite displacement 
(x n ) = (b n + frfi). By considering the system at the lo- 
calized fixed point, A = 0, it is easy to show that - due 
to the logarithmic discretization scheme of NRG - (x n ) 
diverges in the low-energy limit for s < 1 according tc 4 
En • This implies that the average boson numbers 
diverge as well: 



\2„s-l 



(23) 



Importantly, (a z ) in this equation is a temperature- 
dependent quantity which thus depends on n as well, but 
will typically saturate below a temperature T sat . 



A. Truncated mass flow 

Trivially, for Nf, = oo there are no truncation effects, 
and (b n b n ) diverges in the low-energy limit n — > oo. 
For any finite Nj,, the divergence of (b n b n ) is cut-off 
at some n cut , where (b n b n ) reaches a value of order 
O(Nb). In the NRG algorithm, n cu t corresponds to a cut- 
off energy or temperature given by T cut — £n cut , below 
which truncation is relevant. From Eq. (|23p we deduce 
N b oc (a^T^ 1 which yields 



T cut oc N t 



-1/(1-8) 



1(^)1 



2/(1-8) 



(24) 



What happens below T cut ? The sites of the Wilson 
chain with n > n cut no longer behave as free bosonic 
states with energy e„: The low-energy states of the 
shifted oscillators are missing, because the truncated ba- 
sis of the lowest Nb occupation-number eigenstates is too 
small to accommodate them. Qualitatively, the energies 
of the oscillators with n > n cu t are systematically shifted 
upwards. By Kramers-Kronig relations this implies that 
the corresponding real part, Re r„, is modified into Re T„ 
for n > n cut , with Re not reaching Re . 

This is where the mass-flow effect comes into play The 
artificial mass m a (T), given by Re^oo — T n ) from the 
"missing" part of the Wilson chain, will no longer de- 
crease oc T s below T cut , but will instead approach a con- 
stant value m cu t dictated by Re(r oo — T^). This constant 



can be easily estimated: Assuming that the oscillators 
with n > n cut are shifted to infinite energy, it is given 
by m cut = «T c s ut ; corrections to this simple estimate will 
modify the prefactor, but not the power-law dependence 
on T cu t. These considerations imply that the artificial 
mass m a will only flow very little below T cut , i.e., the 
mass flow is truncated at T cut . 

Taking the mass-flow and truncated mass-flow con- 
tributions together, we arrive at an artificial mass 
term in the order-parameter theory, which scales as 
[max(T, T cut )] s . Accounting for different prefactors and 
inserting Eq. (|24[) we arrive at the non-analytic mass 
terms as postulated in Sec. Ill Dl above: 



5 err = / dr4> 2 {r) [i;T s + H(0}r /(1 ~ s) 



(25) 



where w oc N 



-8/(1-8) 



vanishes as Nb — > oo. 



As mentioned in Sec. Ill Dl the non-analytic w term 
asymptotically dominates over the quartic u term for 
s < 1/2 in the mean-field equation of state, Eq. (|19[) . 
such that it leads to the incorrect critical exponents (|2"0)) 
observed in NRG. 

In the power-counting sense, the non-analytic w term 
is marginal for all s: Using Eq. (|12p we have 



dimku] = 1 — di 



= 0. 



(26) 



The system is thus effectively locked at an upper-critical 
dimension. This explains why the NRG exponents (|20|) 
obey hyperscaling, but are nevertheless "trivial" in the 
sense that only simple fractions of s occur. (Note that 
logarithmic corrections do not occur, due to the absence 
of fluctuation effects from w.) 



B. Example: Order-parameter exponent S 

To illustrate the above considerations, wc now discuss 
explicitly the behavior of the magnetization at the critical 
point in a small applied field e as function of temperature 
T, which yields the critical exponent 5 according to e oc 
(j) s as T -> (where = (er z )). 

We discuss the physical behavior of the (j) 4 model © 
in the absence of systematic errors first. The leading 
temperature dependence of the magnetization stems from 
the temperature- induced order-parameter mass m(T). In 
the mean-field regime of s < 1/2, this can be estimated 
from the equation of state: 



= e + m(T)(f> + Zucj? 



(27) 



Since m(T) as T 0, <f> will saturate below T sat , 
with m(T sat ) oc cj) 2 {T = 0). Using m(T) oc T 1 / 2 and 
S = 3, valid in the mean-field regime, we obtain 



e 4 / 3 (s<l/2). 



(28) 
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This argument can be generalized to the non-mean-field 
regime of s > 1/2, using m(T) oc T s , leading to 



= 0.4 



s = 0.6 



t s; 



e 2/(*<5) ( s >l/2). 



(29) 



From Eq. ([27|) one also sees that, above T sa t, </> will follow 
4> oc e/m(T). 

Now we turn to the influence of the systematic NRG 
errors. This requires to discuss the relation between T sa t, 
the physical saturation temperature of cj>, and T cu t, the 
temperature below which Hilbert-space truncation spoils 
the calculation. 

First assume that T sat > T cut . Then the magnetization 



will have taken its T = value 
T cut . From Eq. J23J we deduce 



oc e 1 / 5 already above 



T cu t oc e 



2/[(l- S )«] 



(if T sat > T c 

1 



(30) 



Comparing this with Eq. (|2"9")l . we see that, for small e, 
indeed T sat > T cut consistently holds for s > 1/2. In such 
a case, we expect that the Hilbert-space truncation will 
not significantly influence the T = value of cf>. This is 
in contrast to s < 1/2, where the assumption T sat > T cut 
is found to be violated. 

Now let us discuss the case T sat < T cut , applying to 
s < 1/2, where Hilbert-space truncation will affect the 
impurity magnetization. First, note that the estimate of 
T sa t d28l) needs to be modified due to the mass-flow error: 
m(T) oc T s yields 



e 2/(3s) 



(s < 1/2, mass flow). 



(31) 



At T cut , the magnetization can be estimated as oc 
e/m(T cu t). Again using m{T) oc T s and combining this 
with Eq. (|2"4"|) we have 



-O 1 

+ 

-a 



0.1 



0.01 



-0.0-1 



0.02 



a) J 




r r ' 

> > > > 1 > > />. > 1 > > 


' / 



10 



0.1 













Xn ff. 






/ / I 






\ / f r 






/ 

/ 






/ 

,t , , , 1 , , , , 1 , 






5 10 15 20 



FIG. 2: (Color online) Evolution of the boson occupation 
number, {jr n b n ), and the magnetization, (cr z ), in the spin- 
boson model at criticality in an applied field e, as calculated 
along the NRG run with A = 1, A = 4 and N b = 50. 
a,b) For s — 0.4 (q c = 0.0355465), both quantities saturate 
at the same n, or equivalently, the same temperature T cut 
(dotted vertical lines), illustrating that Hilbert-space trunca- 
tion affects the magnetization. c,d) In contrast, for s = 0.6 
(q c = 0.06702411) the magnetization saturates before T cut 
is reached, and is thus unaffected by truncation. The dark 
dashed lines indicate a divergence (b^bn) oc e^ _1 according 
to Eq. (|23[) ~ this is not realized for s = 0.4 because (cr z } 
is never constant above T cut , but applies to s — 0.6 in the 
interval T cu t < T < T sat where {a z } is constant. 



4>(T cut ) oc 



</>(T cut ) 2s /(i-s) 



(32) 



and 



T cut oc e 2 /( 1+s ) ( s < 1/2, mass flow). (33) 



As argued above, due to the truncated mass flow, m re- 
mains of order T c s ut for all T < T cut . Therefore Eq. ([32] 
implies that 



^^(T cut )« e (1 - S)/(1+S » 



(34) 



Importantly, Fig. [5] illustrates the qualitative difference 
between s < 1/2 and s > 1/2 independent of the knowl- 
edge of the underlying field theory: Hilbert-space trun- 
cation limits the magnetization only for s < 1/2, where 
it consequently spoils the physical critical behavior. 

Upon increasing the applied field e for s < 1/2, T sat 
(f3"Tj) will increase faster than T cut (|33[) . such that the 
correct mean-field behavior is recovered at T = for in- 
termediate values of e. 19 This crossover is fully described 
by the equation of state 



holds down to T = 0, resulting in the critical exponent 
S = (1 + s)/(l — s) as indeed observed. 

The truncated mass flow is shown using NRG data in 
Fig. [2] For s < 1/2 the magnetization <fi = (a z ) indeed 
increases down to Tut, but stops to increase below T cu t 
- this proves both T sat < T cu t and m(T) rj const for 
T < T cu t- The latter reflects the advertised truncated 
mass flow and leads to incorrect scaling. In contrast, for 
s > 1/2 the magnetization saturates already above T cut , 
indicating T sat > T cu t as argued above. 



|e| = t D7V- s/(1 - s V (1+s)/(1 " s) + 3W> 3 (35) 

obtained from minimizing (|19[) for m = 0, i.e., at critical- 
ity, and setting w — wN b s \l — s)/2. A comparison 
of this prediction with numerical data in shown in Fig. [3l 
with excellent agreement except for very small Nb where 
the continuum approximation w.r.t. Nf, is not yet accu- 
rate. 
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FIG. 3: (Color online) Effect of Hilbert-space truncation on 
the evaluation of the critical exponent S in the spin-boson 
model for s — 0.3. The data points have been obtained using 
the variational matrix-product-state approach of Ref. Q11 with 
parameters A = 2, A = 0.1 where a c = 0.0346142; for small 
Nb results from standard NRG are identical up to minimal 
deviations. The lines show the crossover behavior according 
to the equation of state ((35)) with w = 1/50 and 3u — 1/39 
- this gives an excellent description of the data for not too 
small Nb. 



V. DISCUSSION 

We have shown that mass flow, i.e., the fact that 
the order-parameter mass depends on the length n of 
the Wilson chain, is the source of all incorrect expo- 
nents obtained by NRG for the spin-boson model: (i) 
The incorrect finite-temperature exponent x arises be- 
cause the rest of the chain is discarded when calculat- 
ing finite-temperature observables at a temperature T n . 
(ii) The incorrect zero-temperature exponents /3 and S 
arise because the chain is effectively truncated at some 
magnetization-dependent n cut where the diverging oscil- 
lator shifts can no longer be accomodated within the 
finite Hilbert space for each boson site. In the order- 
parameter field theory, (i) and (ii) can be captured as 
artificial non-analytic mass contributions, Eq. (|18[) . 

The fact that the resulting NRG exponents are inter- 
nally consistent and obey hyperscaling can be rational- 
ized by realizing that both errors of bosonic NRG can 
be associated with stationarity along the NRG flow. (A) 
The mass-flow error enforces the temperature-dependent 
renormalized mass to be stationary at criticality (which 
corresponds to naive scaling m a oc T s ). (B) The Hilbert- 
space truncation error cuts the divergence of the boson 



occupation number, making it stationary. Thus, both er- 
rors force the NRG flow to show fixed-point behavior - 
in the NRG sense of stationarity of the many-body en- 
ergies and states along the flow - in a situation where 
the correct physical RG flow does not show stationarity 
in the same sense. This applies both to the localized 
phase of the spin-boson model, where the boson number 
diverges along the Wilson chain, and to the flow towards 
the Gaussian critical point, where the true many-body 
level spacing goes to zero as T — > because of the irrel- 
evance of the quartic interaction^ 

This finally solves the puzzle: NRG produces incorrect 
exponents obeying hyperscaling because its algorithmic 
deficiences force the NRG flow to be stationary, which im- 
plies naive scaling. One concludes that it is difficult with 
traditional NRG to describe critical phenomena which 
violate hyperscaling and are controlled by a dangerously 
irrelevant variable because in such a case the NRG flow 
will not be stationary near the physical RG fixed point 30 
- this has been (partially) illustrated in Ref.l20|. Algorith- 
mic modifications, e.g. the variational matrix-product- 
state approach of Ref. Il3l . can in principle solve the prob- 
lem. 

The present insights show that NRG-based claims of 
violated QCC in other Ising-symmetric impurity mod- 
els with sub-ohmic bosonic bath s 33 ' 34 should be revis- 
ited. Also, implications for the calculation of dynamical 



propertie; 



,35-38 



need to be investigated. In bosonic impu- 



rity models with higher symmetries, QCC was found to 
be fulfilled for a rotor model^ but argued to be violated 
for XY-symmetric spin-boson 1 ^ and SU(2)-symmetric 
Bose-Kondo models. 40 These results suggest that the vio- 
lation of QCC is related to the presence of multiple baths 
(or long-ranged interactions) which couple to different 
non- commuting impurity operators. Clearly, a deeper 
understanding of the relevant non-classical phase tran- 
sitions is desirable. 
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